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Abstract 

We analyze quantum tunneling with the Ohmic dissipation by the non-perturbative 
renormalization group method. We calculate the localization susceptibility to evalu- 
ate the critical dissipation for the quantum-classical transition, and find considerably 
larger critical dissipation compared to the previous semi-classical arguments. 

PACS numbers: ll.10.Hi, 03.65.-w, 73.40.Gk 

1 Introduction 

Caldeira and Leggett proposed a method to derive the quantum dissipative behavior 
from a microscopic theory, which consists of the target system and the environment [1]. 
The action of their model is written as follows, 



S[q,{x a }]=Jdt ^Mq 2 -V (q)+Y: 



(1) 



where q(t) is the variable of the target system in a potential V (q), and x a (t) is the 
harmonic oscillators representing the environment. The target system is coupled linearly 



to each oscillator with strength C a . If we set the parameters m c 



0J r 



[> 0), C a in a 



suitable way, the dissipative term (for the Ohmic dissipation, q) arises in the effective 
classical equation of motion of q(t) after elimination of the environmental variables with 
the proper boundary condition. 
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In this letter, we study the quantum mechanics of this system by the Euclidean path 
integral over q{r) and x a (r). We integrate the variable x a to define the effective action 
for the target system, 

Z = ^fjvqlljvxa = jfjv<l "^ A%1 , (2) 

where S E is the Euclidean action 

Se[q, {x a }} =J dr ^Mq 2 + V (q) + \ m c^a + + ^ C '" X «}> ( 3 ) 

Af and M q are normalization constants and ASfg] is a term generated by the quantum 
effects of environment x a . Caldeira and Leggett studied the influence of AS* on the 
quantum tunneling of the target system q and they found that the quantum tunneling is 
suppressed by the Ohmic dissipation [1]. Following their analysis, many theoretical works 
have been done and the localization (quantum-classical) transition has been suggested to 
occur [2, 3, 4, 5, 6, 7, 8]. The validity of their results have been suggested by elaborate 
experiments [9]. However, the theoretical works so far (instanton, perturbation, etc.) 
depend on some expansions with respect to small parameters which are not always valid. 
For example, for evaluation of the tunneling rate through the double well barrier, 

V (q) = -\Mulq 2 + A g 4 , (4) 

the dilute gas instanton calculation is valid only in A — > region [10]. 

In this letter we adopt the non-perturbative renormalization group (NPRG) method 
[11, 12, 13] to analyze the dissipative quantum tunneling in the Caldeira-Leggett model. 
Since the NPRG is formulated without any series expansion, it is a powerful tool to analyze 
the non-perturbative features of quantum field theory, such as the phase structure of the 
theory [14]. For the quantum tunneling problem, the NPRG method can reproduce the 
exact tunneling rate through the double well barrier (Eq. (4)) in the wide parameter Ao 
region [10, 15, 16]. Therefore the NPRG is expected to be effective in the non-perturbative 
analysis of the dissipative quantum tunneling. We can discuss the possibility of the 
localization (quantum-classical) transition by investigating the "phase structure" of the 
dissipative quantum mechanics. The NPRG also helps us to interpret the effects of AS 
as effective infrared or ultraviolet cutoffs for the quantum fluctuation, and we can readily 
understand the dissipation effects as suppression or enhancement of the quantum features 
of the effective target system [17]. 

2 Dissipative effective action 

The effective interaction AS* defined in Eq. (2) takes the following form, 

A%] = - J drj ds q(r) a(r - s) q(s). (5) 

The non-local coupling coefficient is given by 

a ( T - s) = f—Y* -3 \ e Mr-s) 

{ } J 2tt ^ 2m a ^ + ool 

7^ M e-^', (6) 
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where we have introduced the spectral density function J{uj) characterizing the environ- 
ment, 



J(u) 



E 




(2ir) 5{u) — uj a ). 



(7) 



Am a uj a 



If we set J{uj) = i] u, the Ohmic dissipation term — r\q arises in the effective classical 
equation of motion of q(t) with the proper boundary condition of the environmental 
variables [5, 7]. 

Generally AS consists of the local component ASl and the non-local component ASnl- 
We identify the component AS N l as the dissipation term, since it corresponds to the odd 
power term in the Fourier transform and then it is related to the time reversal symmetry 
breaking. On the other hand, the component ASl consists of the even power term in the 
Fourier transform and then it does not contribute the time reversal symmetry breaking 
in the effective classical equation of motion, that is, it is irrelevant to the dissipation. 
Therefore, we introduce a suitable conterterm in the Euclidean action S E to cancel the 
component AS L [1]. 

For the Ohmic dissipation J{uj) = r/u, we identify the dissipative part as follows: 



It should be noted that the dissipation term, ASnl[?], never breaks the time reversal 
symmetry in the action level. The actual dissipative effect (—rjq) arises only in the clas- 
sical equation of motion level [5, 7]. Using the Fourier transform of q(r), AS'nl is also 
represented as [7] 



Note that rj is a dimensionful parameter, [77] = pVLj] = [mass] x [time] 1 . We study the influ- 
ence of the dissipation term ASnl on the quantum behaviors of qij). 

3 NPRG equation with the Ohmic dissipation 

Now let us proceed to the NPRG analysis of the Caldeira-Leggett model [18]. Origi- 
nally the NPRG method has been formulated and used mainly in the statistical mechanics 
or the quantum field theory to analyze the critical phenomena [11, 12, 13, 14]. Recently, 
the NPRG method has been found even effective in the quantum mechanics particularly 
for the non-perturbative analysis [10, 15, 16]. 

We derive the NPRG equation for the quantum mechanics with the dissipation term. 
In the NPRG method, the theory is defined as an effective theory with a high frequency 
cutoff A, and is described by the Wilsonian effective action Sa[<?], which changes with 
respect to the scale A. We now employ the local potential approximation (LPA), 






(9) 



SM = J 



ds 



(q(r)-q(s) f 



+ V K {q) 



(10) 



T — S\ 2 
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where the effective action is limited to have only local potential term, the Wilsonian 
effective potential V\(q), in addition to the fixed kinetic term and the fixed dissipation 
term. Both the particle mass M and the dissipation strength 77 are constant and only the 
Wilsonian effective potential V\(q) is changed as the scale A is lowered. Note that the 
LPA is the leading order of the derivative expansion of S\[q] and has nothing to do with 
the prescription that we ignore the local part ASl in the component AS. The LPA means 
that any quantum correction to the derivative coupling is ignored, but does not mean that 
any quantum correction from the derivative coupling is ignored. We carry out the path 
integration over the highest frequency degrees of freedom of q(u), A — AA < < A, up 
to the one-loop order, 

V A - AA (q) = V A (q) + h / a A aa ^ log (Mu 2 + V \u\ + ^) , (11) 

where we have omitted a proper constant term which should make the argument of the 
logarithm to be dimensionless, since it contributes only to the g-independent part of the 
effective action. Taking the limit AA — > in Eq. (11), we have the NPRG equation 

dV\ h . , / 1 77 11 d 2 V A \ 

A W = -2? a1o H 1 + am + a5mT^J- < 12 > 

which should be called the LPA Wegner-Houghton equation with dissipation. This is 
because in the 77 — > limit Eq. (12) becomes the ordinary LPA Wegner-Houghton equation 
[12, 13], which describes the change of the Wilsonian effective potential V A (q) with respect 
to the scale A. The extra term rj/(AM) originates from the dissipation term ASnl- The 
NPRG equation (Eq. (12)) is obtained with the approximation LPA; the leading order 
of the derivative expansion of S\ [q] ■ The derivative expansion does not rely on any small 
parameter which controls the series expansion, and therefore the effectiveness of the NPRG 
equation is expected to be free from the smallness of the parameters in the system. 

4 Analysis of Dissipative Quantum Tunneling 

We regard the system variable q as a macroscopic collective coordinate and discuss 
how the quantum mechanical behavior of the macroscopic coordinate q is affected by the 
environmental effects. As a typical phenomenon, the quantum tunneling of q has been 
analyzed by many authors [1, 2, 3, 5]. We set the bare potential of q as the double well 
type (Eq. (4)). In this system with vanishing 77, the coordinate q tunnels through the 
potential barrier and it oscillates between two wells. The main question is whether it also 
oscillates or not even in the 77 7^ quantum dissipative dynamics. Such an oscillating 
behavior is particularly called the macroscopic quantum coherence and the inclination of 
ceasing the oscillation corresponds to decoherence [6]. 

Now we briefly summarize the previous results obtained for the macroscopic quantum 
coherence in double well potential systems with the Ohmic dissipation. Usually the first 
energy gap AE = E\ — Eq is calculated as the physical quantity because it corresponds 
to the tunneling amplitude between two wells. Caldeira-Leggett evaluated it by using the 
semi-classical (instanton) approximation which is valid for A — > and the perturbation 
with respect to 77, and then they found the dissipation term ASnl suppresses the quantum 
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tunneling [1]. Renormalization group analyses have been done for the instanton gas 
system within the dilute gas approximation which is valid also for A — > and 77 — ^ 
region [2, 3]. They predict a remarkable phenomenon that AE vanishes at a critical value 
f] c = 2kTi\q/ (MuJq), where the decoherence occurs. This phenomenon is often called the 
quantum- classical transition because the system looks classical in a limited sense that the 
quantum tunneling does not occur [19]. This is nothing but the spontaneous Z% symmetry 
breaking. In these works the energy gap AE was treated as the order parameter of the 
transition. 

These results are obtained by using the semi-classical approximation and/or the per- 
turbation theory with respect to 77. Their reliability depend on the smallness of the 
couplings Ao and rj. Therefore, to get more general and reliable results free of such limi- 
tation, we must employ an analyzing tool which does not need the series expansion with 
respect to any couplings. We expect our method using the LPA Wegner-Houghton equa- 
tion with dissipation (Eq. (12)) will work best for the large coupling region, because the 
employed approximation, LPA, is the leading order of the derivative expansion and does 
not depend any small parameter. In fact, it has been found that the NPRG analysis of 
quantum tunneling based on the LPA Wegner-Houghton equation works very well in the 
wide parameter (A ) region [10]. 

We analyze the dissipative quantum tunneling by solving the NPRG equation (Eq. 
(12)) numerically. This equation is a two-dimensional partial differential equation for the 
Wilsonian effective potential V\(q) with respect to A and q. The initial condition of the 
potential is the bare potential Va. (?) = — \Muj^q 2 + A g 4 at the initial cutoff A . We 
solve the differential equation toward the infrared limit A — > and finally obtain the 
physical effective potential V cS (q) = ^Mu 2 s q 2 + A eff g 4 + • • •. We can exploit the physical 
information of the quantum system from the effective potential V e g(q). First of all, the 
frequency squared starting with the negative value u\ Q = — uJq finally reaches a positive 
value (for small rj), u 2 s > 0, that is, the effective potential V e s(q) becomes a single well 
form. It means that due to the quantum tunneling, the Z 2 symmetry (parity) does not 
break spontaneously, and the expectation value of q is vanishing. The tunneling effects 
are automatically incorporated when we integrate (solve) the NPRG equation toward the 
infrared. 

It should be noted that for non-vanishing r], the simple equivalence between the energy 
gap AE and the effective frequency u; cff , AE = hu cS [10], does not hold. The long range 
(r — > 00) behavior of the two point function 

lim (Q\Tq(r)q(0) \Q) = lim / ^e wi "-— ^ — (13) 

is actually independent of uj c k and is determined by rj with power damp behavior even in 
case of infinitesimal r]. This singular behavior comes from the nature of the environmental 
degrees of freedom where the harmonic oscillator frequencies are assumed to distribute 
continuously up to zero (cu a = 0). Even if there is an infrared cutoff, it does not change the 
situation much and the lowest frequency environment dominates the long range correlator 
of the target system. In this sense the semi-classical arguments calculating the first 
energy gap AE = E\ — Eq as the order parameter of the quantum-classical transition in 
the quantum dissipative systems seems somewhat doubtful, because the correspondence 
between the tunneling amplitude and the first energy gap AE comes form the long range 
behavior of the transition amplitude. 
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Here, instead of analyzing the long range correlator, we adopt another physical quan- 
tity to describe the possible quantum-classical transition, the localization susceptibility 
We add a source term 



,T/2 

Jq T (0) = J dr q{r 

J-T/2 



(14) 



to the Euclidean action Se (Eq. (3)). Then the localization susceptibility x is defined by 
d(q(r))j 



X = 



dJ 



j=o 



lim - d &(°» J 



d,J 



j=o 



,. 1 d 2 \ogZ 
lim — 

T^oo T dJ 2 



(15) 



J=0 



and it is exactly given by the effective potential as follows: 



X 



d 2 V eS (q) 



<?=9n 



(16) 



where g min is a minimum of V e g(q) and it is actually vanishing in the sub-critical region. 
These definitions are valid also for non-vanishing 77. If we find a divergent and scaling 
behavior of x toward rj = rj c , we may conclude that the localization transition occurs at 
rj c and it is the second order phase transition. 




Figure 1: The effective frequency u e g with Figure 2: Critical scaling fit of localization 
initial cutoff A = lC^f/kuo]- susceptibility x f° r -^o — 1.0[M 2 a>o/fr]. 

The results are shown in Fig.l for several values of Ao- In previous works for rj = 
system [10, 15, 16], the NPRG analyses have been found to work excellently in these Ao 
region (Ao > 0.1[M 2 ci;q/^]) while the dilute gas instanton approximation does not work 
at all there. Therefore, we expect the NPRG analysis also works well for these Ao values 
even in 77 7^ systems. We find that for every value of A , uj ef[ decreases as 77 becomes 
large. This r\ and Ao dependence of uj c q is naturally understandable by considering the 
dissipation effect in Eq. (12). The quantum fluctuation of the dynamical variable q 
is suppressed bellow the effective infrared cutoff c^ir = rj/M. Therefore, the larger 77 
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causes the stronger cutoff effect and results in the smaller oj e g. For larger Ao, the effective 
frequency of the system is larger, and therefore larger r\ is required to effectively cutoff 
the quantum fluctuation. These behaviors are qualitatively consistent with those of AE 
obtained by the instanton approximation. 

Then, what happens for large rj case? The expected phenomenon is the quantum- 
classical transition characterized by complete disappearance of the tunneling. The NPRG 
method may not work precisely in the u cS /uj — > region because Eq. (12) becomes 
singular there and the numerical error increases. Therefore, according to the standard 
technique of analyzing the critical phenomena, we evaluate the critical dissipation i] c from 
the diverging behavior of the localization susceptibility. We fit x(v) with the critical 
exponent form 

x = C\r l -r lc \-\ (17) 

We show an example of fit for Ao = 1.Q[M 2 oJq/K\ in Fig. 2. We conclude that the localiza- 
tion susceptibility exhibits a divergent behavior with a power scaling, which indicates the 
second order phase transition of localization. We obtain the critical dissipation r] c and the 
critical exponent 7 for these Ao, which are plotted in Fig. 3. The previous analysis using 
the dilute gas instanton gives a simple relation i] c = 27c?iXo/(MuJq) which is also shown 
in Fig. 3 for comparison. The NPRG results for r\ c are systematically larger compared to 
those of the instanton. As for the critical exponent 7, we observe the universality property 
except for the smallest A = O.lfM 2 ^/^,], where NPRG results may not be reliable. 
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Figure 3: A dependence of critical dissipation rj c and critical exponent 7. 



The NPRG works very well in the large Ao region [10, 15, 16]. Strictly speaking, 
tunneling phenomena are no longer characteristic in such parameter region. It can be 
easily seen by evaluating a dimensionless quantity r = 8\^2h\o / (M 2 u>q) , which is the ratio 
of the zero-point energy of the particle at each well with 7] = to the height of the inter- 
well barrier. Tunneling phenomena are dominant for r < 1 case. For Aq = 0.1 [M 2 Uq/H], 
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r ~ 1.1 is obtained. Therefore, our analysis in the large Ao region (Ao > 0.1[M 2 u1q/?i]) 
rather claims suppression of the quantum oscillation over the small inter-well barrier. The 
critical phenomenon that we observed here is induced by a destruction of the quantum 
coherence and can be interpreted as spontaneous Z 2 symmetry breaking-like phenomenon 
in quantum mechanics. 

5 Summary 

The Caldeira-Leggett model for dissipative quantum mechanics was analyzed by 
means of the NPRG method. We derived the LPA Wegner-Houghton equation with 
dissipation, which is the local potential approximated NPRG equation with the Ohmic 
dissipation term. We applied it to analysis of dissipative quantum tunneling (or dis- 
sipative quantum oscillation) and investigated whether the quantum-classical transition 
occurs or not. Since the first energy gap AE = E\ — E is unsuitable as the order pa- 
rameter of the dissipative phase transition, we calculated the localization susceptibility x 
from the effective potential V cS (q). The observed x diverges toward rj = r] c (\ ) with the 
critical scaling, and therefore we concluded that ?7 c (Ao) is the critical dissipation for the 
quantum-classical transition and this transition seems the second order phase transition. 
The values of the critical dissipation 77 C (A ) are rather larger than the values obtained by 
the semi-classical analyses, where the first energy gap AE is used as the order parameter 
of the quantum-classical transition. As for the non-Ohmic dissipations, the analogous 
systematic study using the NPRG equation will be reported elsewhere [17]. There must 
remain, however, many fundamental and subtle problems to be studied about the notion 
of quantum-classical transition and the critical dissipation. The NPRG method will be 
effective also for those problems. 
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